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Abstract 

Cln The axisymmetric kinematic dynamo problem is reconsidered and a number of 

open questions are answered. Apart from axisymmetry and smoothness of data and 
solution we deal with this problem under quite general conditions, i.e. we assume a 
compressible fluid of variable (in space and time) conductivity moving in an arbitrary 
(axisymmetric) domain. We prove unconditional, pointwise and exponential decay of 
magnetic held and electric current to zero. The decay rate of the external (merid- 
ional) magnetic field can become very small (compared to free decay) for special flow 
fields and large magnetic Reynolds numbers. We give an example of that. On the 
other hand, we show for incompressible fluids of constant conductivity in balls that 
the meridional and azimuthal decay rates coincide with the well-known poloidal and 
toroidal free decay rates, respectively. 

Key Words: Magnetohydrodynamics, dynamo theory, axisymmetric theorem. 

1 Introduction 

The axisymmetric (or Cowling's) theorem is the oldest and most renowned of all antidy- 
namo theorems which generally preclude under certain assumptions the maintenance of 
the magnetic field against ohmic loss by the dynamo process. Here, the assumption is 
axisymmetry (of magnetic field, flow field, conductivity distribution and shape of the con- 
ductor). This theorem had considerable significance in the early days of dynamo theory, 
where convincing analytical arguments, computational or even experimental evidence for 
the dynamo process were yet far away. In view of nearly axisymmetric magnetic fields 
of the earth and the sun this theorem called in question the dynamo process as such. 
No wonder that after its very first provisional formulation by Cowling (1934) researchers 
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in the field tried to improve the theorem with respect to greater mathematical rigor, 
stronger statements or weaker assumptions. The history of this endeavour has already 
been recorded in sufficient detail elsewhere. I mention only the excellent articles by Ivers 
& James (1984), who cover the situation up to 1983, by Fearn et al. (1988), who note 
in particular that this endeavour is not yet finished, and by Nunez (1996), whose review 
presents the state of the art and emphasizes especially the role of Lortz and his coworkers. 
So, I refrain from giving another detailed account of the history but distinguish roughly 
three phases in the development of the axisymmetric theorem. 

In the first phase the stationary problem with solenoidal flow field and constant con- 
ductivity was in the focus. Cowling himself demonstrated the impossibility of a purely 
meridional magnetic field sustained by a purely meridional flow field by means of his now 
famous 'neutral point argument'. Backus & Chandrasekhar (1956) generalized this result 
by admitting meridional as well as azimuthal components for the magnetic field and the 
flow field, and they replaced the somewhat informal neutral point argument by mathe- 
matically more rigorous maximum principles for elliptic equations. Still later Lortz (1968) 
realized that these elliptic equations are incompatible with nontrivial solutions even for 
non-solenoidal flow fields and variable (in space) conductivity. 

In the second phase the more realistic time- dependent problem came into focus but 
still under the proviso of a solenoidal flow field and constant conductivity. With these 
assumptions energy methods and variational inequalities can be applied to prove expo- 
nential decay in the energy norm of the meridional scalar and, in the absence of the 
meridional field, of the azimuthal field. Backus (1957) and Braginskii (1965) did so and 
gave, moreover, for balls explicit (albeit different) decay rate bounds on these quantities. 
Note that Backus (1957) considered already at that time the meridional problem without 
the solenoidal-flow and constant-conductivity asumptions and found exponential decay; 
however, in order to evaluate the variational problem, he had to introduce unsatisfactory 
ad-hoc assumptions. 

The third phase was initiated by a paper of Todoeschuck h Rochester (1980), who, in 
the view of Saturn's almost axisymmetric magnetic field, argued that Cowling's theorem 
might have loopholes, in particular, that non-stationary axisymmetric magnetic fields 
could be generated by non-solenoidal flows in a compressible fluid. Several researchers, 
however, dismissed this possibility and tried to close the loopholes; so, the third phase was 
characterized by attempts to prove the axisymmetric theorem for non-solenoidal flows and 
variable conductivity. The first such attempt was due to Hide & and Palmer (1982), who 
demonstrated monotonous decay of the meridional scalar in the maximum norm. Their 
arguments were ingenious but, from a mathematical point of view, not rigorous (cf. Nunez 
1996). Lortz & Meyer-Spasche (1982a) gave a rigorous proof of this result using a clever 
combination of parabolic and elliptic maximum principles. Ivers & James (1984) combining 
these maximum principles with skillfully chosen comparison functions strengthened this 
result in that they showed exponential decay to zero together with explicit decay rate 
bounds. 

Less is known about the azimuthal field. First of all the azimuthal problem so far has 
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exclusively been considered in the case that the meridional field has already died out and 
does not provide a source term for the azimuthal field. With this assumption Lortz & 
Meyer-Spasche (1982b) proved the azimuthal field to be bounded in an integral norm for 
all times in terms of its initial value. Their method of proof depends of the existence of 
positive solutions of an auxiliary problem, which requires some mathematical effort. Ivers 
& James (1984) obtained the same result using some more intuitive arguments. Finally, 
using heavy mathematical machinery, Lortz, Meyer-Spasche & Stredulinsky (1984) proved 
in a plane model problem exponential decay of the azimuthal field to zero even in the 
maximum norm. 

Despite all these achievements there are still a number of shortcomings that need to 
be removed before the axisymmetric theorem can be regarded as fully established: 

(i) In an axisymmetric setting the evolution of the meridional field decouples from 
that of the azimuthal field (which is the essence of the axisymmetric theorem), but not 
vice versa. Nevertheless, the evolution of the azimuthal field so far has exclusively been 
considered without meridional source term. So, a proof of the axisymmetric theorem using 
the full (coupled) set of dynamo equations is still missing. 

(ii) Meridional decay results refer so far to the meridional scalar, for which a tractable 
equation is at hand, and not to the meridional field itself, with two notable exceptions: 
Hide (1981) argued that a signed version of the total magnetic flux leaving/entering the 
conductor has to decay in the case of axisymmetry; his arguments, however, are not 
rigorous. Ivers & James (1984) appeal to Schauder estimates for elliptic and parabolic 
equations to relate the decay of the meridional field to that of the meridional scalar; these 
estimates, however, are not uniform in space; they fail, in particular, at the magnetic axis 
and at the conductor's boundary leaving thus open - however unlikely - the possibility of 
magnetic regeneration there. 

(iii) Even in the absence of the meridional source term, unconditional decay to zero 
of the azimuthal field is not yet established. Ivers & James (1984) discuss a 'scenario for 
decay to zero' based on additional ad-hoc assumptions, whereas the result of Lortz et al. 
(1984), although highly suggestive, does not make a statement about the axisymmetric 
problem. 

(iv) The only explicit decay rate bounds valid for general flow fields and conductivity 
distributions refer to the meridional field and are due to Ivers & James (1984). These 
bounds become extremely small already for moderate values of the magnetic Reynolds 
number R m . For instance, considering the earth, one finds for R m = 10 2 (which is a rough 
estimate) magnetic decay rate bounds such that the corresponding decay times greatly 
exceed the lifetime of the earth. So, the question arises how close these bounds are to 
actual decay rates. 

(v) In the case of solenoidal flow and constant conductivity the optimal decay rates 
are those of axisymmetric free decay. In balls these can be compared to the well-known 
poloidal/toroidal free decay rates. Calculations by Backus (1957) suggest here a difference 
whereas Braginskii (1965) notes (without proof) no such difference. 

We prove in the following unconditional, pointwise and exponential decay to zero of 
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the meridional and the azimuthal fields satisfying the coupled dynamo equations, thus 
settling, finally, problems (i) to (iii) of the above list. Problem (iv) is resolved by example. 
We prove for a special piecewise constant flow field that the meridional decay rate shrinks 
exponentially fast to zero with respect to R m . This behaviour does not occur in simpler 
'one-dimensional' flows. For purely radial or purely non-radial flows it is proved that 
meridional decay rates do not drop significantly below the corresponding free decay rates, 
irrespective of R m . Concerning the last problem (v) we derive the governing equations of 
meridional and azimuthal free decay. In balls we find the corresponding decay rates equal 
to those of poloidal and toroidal free decay. We thus verify the assertions of Braginskii, 
but make also a comment on Backus' differing result. 

In proving these results we make use of all the mathematical tools already applied by 
previous authors, in particular, we use variational methods as in (Backus 1957), compar- 
ison functions as in (Ivers & James 1984) and positive auxiliary solutions as in (Lortz et 
al. 1984). The main idea, however, which allows to go beyond these authors, is a consis- 
tent use of the 5-dimensional formulation of the axisymmetric problem. This formulation 
eliminates the coordinate singularity inherent to the 3-dimensional axisymmetric formu- 
lation and makes the problem accessible to more or less standard techniques for partial 
differential equations. For instance, Galerkin representations of solutions are essential in 
the present context. These are standard for the azimuthal field and have recently been 
provided for the meridional scalar (Kaiser & Uecker 2009). For solutions in this form we 
can use a 'higher-order-decay theorem' relating the decay of the solution to the decay of 
its higher derivatives (Kaiser 2012b). This result allows, in particular, to relate the decay 
of the meridional scalar to that of the meridional field, and, moreover, to relate integral 
decay to pointwise decay. 

Concerning the organization of the material I have tried to make the paper accessible 
also to those readers who are less inclined to mathematics but are interested in the status 
of the axisymmetric theorem. So, some more formal aspects of the paper, some proofs and 
some auxiliary calculations have either been omitted (and published elsewhere) or shifted 
to appendices, and the trustful reader might concentrate on the conclusions in the main 
text. In more detail the paper is organized as follows: section 2 presents the formulation 
of the axisymmetric dynamo problem in 2, 3 and 5 dimensions, which emphasize different 
aspects of the problem. In section 3 the higher-order-decay theorem is presented and 
then applied to the meridional decay problem whereas in section 4 results by Lortz et 
al. (1984) are adjusted to the axisymmetric situation and then applied to the azimuthal 
decay problem. A more technical part of this section has been shifted to appendix A. 
Section 5 is devoted to the construction of meridional sub- and supersolutions to obtain 
explicit upper and lower decay rate bounds on the external magnetic field for certain flow 
fields and section 6, finally, presents sharp meridional and azimuthal decay rate bounds for 
incompressible fluids of constant conductivity. Technical details have again been shifted 
to the appendices B and C. 
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2 Formulation of the axisymmetric dynamo problem 



In the framework of magnetohydrodynamics the kinematic dynamo problem is the follow- 
ing initial value problem in all space (Backus 1958, Moffatt 1978) 



d t B 



-V x (r/V x B) + V x (v x B) , V B 
VxB = 0, V • B = 
B continuous 
B(x,-) -»• 
B(-,0)=B , V-B = 



in G x R+, 
in G x R + , 
in R 3 x R+, 
for |x| — > oo, 
on G x {t = 0}. 



(2-1) 



Here, the induction equation (|2.ip i describes the generation of the magnetic field B by 
the motion (with prescribed flow field v) of a conducting fluid (with magnetic diffusivity 
rj > 0) in a bounded region G C H 3 . Outside the fluid region there are no further sources 
of magnetic field. Thus, B matches continuously to some vacuum field in G := H 3 \ G 
that vanishes at spatial infinity. 

The central simplifying assumption of the present paper is axisymmetry of all variables 
appearing in (|2,ip including the shape dG of the conductor. Using cylindrical coordinates 
(p, <j>, z) with e z pointing in the direction of the symmetry axis S and (p, 4>) being po- 
lar coordinates in the planes perpendicular to S, this assumption implies the following 
representation of the solenoidal field B 



B = VM x + AVcj) 



d z Me p + 



d p Me, 



+ -ie, = B m + B fl . 



P 



(2.2) 



The meridional field B m is confined to planes <fi = const and thus everywhere perpendicular 
to the azimuthal field B a . The corresponding scalars M and A depend on p, z and t but 
not on (j>. 

Inserting (|2.2|) into (|2.ip i„ and using the abbreviation E := rjV x B — v x B yields 
comp onentwise 

d g (dtM + pEj,) = d^E,, 
d t A + 8 z E p - d p E z = 0, 

d p {d t M + pE^) = d^Ep. 

The right-hand sides of (|2.3p i 3 vanish by assumption and one concludes 



(2.3) 



d t M + p(nV x (VM x + AX7(p) - v x (VM x V^ + AVcp)), = c(t), 



(2.4) 



where c(t) is an undetermined function of t. Evaluating the left-hand side of (|2.4p and 
replacing M — J 1 c(r) dr by M, which does not affect the meridional field, we arrive at 



dtM - ??A^M + v p dpM + v z d z M = 0, 



(2.5) 
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where A* denotes the elliptic operator d 2 — 1/p d p + d z . 

On the other hand, evaluating (|2.3p ? yields an evolution equation for the azimuthal 
scalar A with a source term depending on M: 

d t A - P d p {^ d p A) - d z (r, d z A) + p d p (^A)+d z (v z d z A) = - p d p (^-)d z M + d zV<t> d p M. 

9 '' (2.6) 

Similarly but simpler, inserting (|2.2p into (|2.ip 9 n yields in the vacuum region 

A*M = (2.7) 

and d p A = d z A = 0, which implies by the asymptotic condition (|2.ip ,i. 

A = 0. (2.8) 

The matching condition (|2.ip ^ implies for M the matching condition 

M, d p M, d z M continuous (2.9) 

and for A by (|2.8p the boundary condition 

A\ dG = 0. (2.10) 

The asymptotic condition (|2.1P a reads for M 

-d p M(p,z,-) -»• 0, -8 z M(p,z,') for p 2 + z 2 -»• oo, (2.11) 
P P 

where this notation implies that the convergence is uniform on compact sets of the variables 
represented by dots. The coordinate singularity at p = inherent to the representation 
(|2.2p requires additional conditions on M and A to avoid spurious solutions: 

M(p, ■, •) -> , d p M(p, ■, •) = O(p) , 8 z M(p, •, •) = 0{p) for p -> (2.12) 

and 

A(p,-,-) = 0(p 2 ) for p^0. (2.13) 

The conditions (|2.12p ^ / . ensure a finite magnetic field on the symmetry axis S and im- 
ply, moreover, the finite limit linip_>o M(p, •, t) =: Ms(t). As Mg(t) does not affect the 
meridional field we set it to zero. Condition (|2.13p ensures a differentiable (with respect 
to p) magnetic field component that vanishes at S. Finally, H 3 is replaced by the 
half-plane H := {(p,z) £ K + x M} and G and G by the 'cross-sections' G2 := {(p, -z) £ 
H : x(p, 0, 2) £ G} and Gi := H \ G2, respectively. 

So, summarizing the foregoing results, a solution of the axisymmetric dynamo problem 
consists in a couple (M, A) with M satisfying (|2.5p in G2 x H + with initial value Mq , (|2.T[) 
in G2 x R + and (12. 9p . (12. lip and (|2.12|) in i? x R + and, furthermore, with A satisfying 
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(|2.6p and (|2.10p in G2 x H+ with initial value Aq. The initial values themselves have to 
satisfy (|2.12p and (|2.13p and typically some 'compatibility conditions' (cf. Kaiser 2012a). 

This formulation in terms of the potentials M and A and the spatial variables p and 
z may be called minimal since it completely resolves the condition of axisymmetry. Most 
important, it yields an evolution equation for M decoupled from A and without zeroth- 
order term. However, there are also drawbacks. The divergence character of the azimuthal 
equation is veiled and, more important, we have to deal with singular coefficients in both 
equations. In fact, more redundant formulations turn out to be useful using cartesian 
coordinates in IR 3 and even R 5 and using modified potentials. 

The 3-dimensional formulation is straightforward. With x = (x, y, z) G H 3 and the 
abbreviation p := (x 2 + y 2 ) 1 ^ 2 the meridional problem can be summarized as follows0 



d t M 



Vp 

r/AM + 2r] — ■ VM + v • VM 
P 

Vp 



P 

M and VM continuous 



, VM(x,y,- 



P 



(2.14) 



in G x R+, 

AM: - 2— • VM - in G x R+, 

in M 3 x M + , 
0{p) for p -»• 0, 

VM(x,')->0 for |x| -)■ 00, 

M(-,0) = M onGx{t = 0}. 

and in (|2.14p i ? we made use of 



Here, V denotes the Cartesian gradient vector (d x ,d y ,d z 
the identity for axisymmetric functions, 



9 2 M 



1 



AM 



Vp 
2 — • VM. 



P P 

Concerning the azimuthal problem we introduce the variable A := A/p 2 and use the 
identity for axisymmetric functions 



to obtain 
d t A 



V- 



l -d p (^d p {p 2 A 



j,V(p 2 A 



d z {7 1 d z A)=V 



+ V • {vA) = V • ( B r 
A = 
^(•,0) =A 



P 



±V(p 2 A) 



in G x R + , 
on x M + , 
on G x {t = 0}. 



(2.15) 



Note that a solution of (12. 15H that is well-defined on the symmetry axis automatically 
satisfies (|2.13p . 



lr To keep the notation simple we do not use different symbols for the same function depending on 
different coordinates. The arguments should be clear from the context. 
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The 3-dimensional formulation of the azimuthal problem makes the divergence char- 
acter manifest; the coefficients, however, are still singular. This problem can be cured 
by a five-dimensional formulation. Let x = (x\, X2, £3, 24, £5) G K 5 , V the corresponding 
cartesian gradient and A the 5-dimensional Laplacian. Identifying p 2 with Yli=i x i an d 
z with £5 and defining G5 := {x 6 M 5 : (p, z) € G2}, axisymmetric functions in H 3 can 
be considered as axisymmetric functions in H 5 . Introducing the variable Ai := M/p 2 and 



the 5-dimensional meridional flow field v 
axisymmetric identity 



(5) ._ 



Y,i=i( v p/P) x i e i + ^ e 5> an d observing the 



A< f M=-d p {p'n i , 



M 



+ d z z M 



the meridional problem takes the form 
d t M 



r)AM + v$ ■ VM + 2 — M = 
AX = 



in G5 x R + , 
in G5 x R + , 
in IR 5 x 1R + , 
for |x| — > 00, 
on G 5 x {t = 0}. 



(2.16) 



.M and VAi continuous 
M(x,-) -> 
M(-,0)=M 

Note that for a smooth axisymmetric flow field v p /p is well-defined on S; thus, no singular 
coefficients appear any more in (|2.16p . No axis-condition on Ai is anymore necessary; 
condition (|2.12j) is automatically satisfied for any well-defined solution of (|2.16p . Moreover, 
outside the conductor, Ai is a harmonic potential (in 5 dimensions); using this information 
(cf. Kaiser & Uecker 2009, Appendix C) condition (I2.16h /i implies, in particular, 

-4h 



V7W(x, 



0(|x| 



for Ixl 



00, 



which again implies (12. lip . Note, however, that for these advantages we had to pay a 
prize, viz. a zeroth-order term in ()2.16j) i . 

Concerning the azimuthal problem we make use of the axisymmetric identities 



pd p (^d p A) +d z {r,d z A) 



d p r) 



and 



to obtain 



pd p (^A) + d z v z A = p 2 V.(v% 4) "2-^ 



d t A - V ■ (rjVA) + V • (v£U) - 2 

= dvf -VM + 2 8 * v + 
A = 
A(-,0) = Ao 



d p r] 



+ 



A 



P 



Ai 



in G 5 x R + , 
on dG§ x R + , 
on G 5 x {t = 0}, J 



(2.17) 



S 



where we used the abbreviation dvi 5 '* := Ylt=o ®z v <l> i x i/p) e i ~ P^p{ v 4>l 'p) e 5- On the 
premises of smooth (up to second order) axisymmetric data r/ and v, all coefficients in 
(|2.17p are now well-defined on S (cf. Ivers & James 1984). But, again, we have bartered for 
this advantage a zeroth-order term that destroys the divergence character of the left-hand 
side in (|2T7|l i. 

As far as the various forms of meridional and azimuthal evolution equations have 
been considered by previous authors, these coincide, of course, with ours. Concerning the 
asymptotic conditions at S and at infinity, however, there are some minor differences. For 
example, Ivers & James (1984) use the condition 

M(p,z,-) =0((p 2 + z 2 y 1/2 ) forp 2 + z 2 ^oo (2.18) 

instead of ()2.11|) . Stredulinski et al. (1986) called condition (|2. 18[) in question. The 
discrepancy, however, seems to be due to a confusion about the axisymmetric and the 
plane-symmetric problem (cf. Kaiser & Uecker 2009, Remark B.l). In fact, in view of the 
five-dimensional formulation conditions fj2 . 1 1 f) and fj2. 18f) are equivalent for solutions of 
()2.7p . A similar remark holds for the axis-condition (|2.12p . which is typically replaced by 
previous authors (as far as they specify at all an axis-condition) by the slightly weaker 
condition 

M(p,-,-) = 0{p 2 ) for p^O. (2.19) 

In fact, ()2.19p is enough to establish the equivalence between the two- and the five- 
dimensional problem and hence the equivalence between (I2.12P and (I2.19P for meridional 
solutions. 

Concerning the flow field some authors require the normal component to vanish on 
dG, others do not. In fact, some results depend on such a condition, others do not. In the 
following we will specify this condition where needed. 

A final remark concerns the question, to which quantities axisymmetry has exactly to 
apply. In view of eqs. (|2.3p Lortz (1968) required axisymmetry of the electrical quantities 
B and E but not of v or r]. Todoeschuck & Rochester (1980) critized this approach as a 
disguise of the usual requirements. Lortz's assumptions are indeed enough to derive eqs. 
(I2.5P and (I2.6P even for non-axisymmetric coefficients v and r/. However, evaluating the 
conditions d^Ep = d§E z = on their part yield the additional equations 

d^T] d p A - d^Vp A - d z M = 0, 
d^rj d z A - d<f,v z A + d^v^ d p M = 0. 

These equations were trivially solved by axisymmetric coefficients (and this solution To- 
doeschuck & Rochester probably had in mind). Otherwise, for prescribed non-axisymmetric 
coefficients, (|2.20p represent additional constraints. The combined system (|2.5p . (|2.6p and 
(I2.20P for axisymmetric scalars M and A, however, does not seem to be well-posed; in 
particular, general solvability cannot be expected. So, the present treatment of the ax- 
isymmetric problem cannot dispense with the assumption of axisymmetric coefficients. 




9 



3 Decay of meridional field and azimuthal current 



In this section the exponential decay of the meridional scalar M in the maximum norm 
is related to exponential decay with (almost) the same decay rate of the meridional field 
B m and of the azimuthal current J a . As discussed in the introduction pointwise decay of 
M, for instance in the form 

\M(p,;t)\<Cg(p)e~ dot , p>0,t>0 (3.1) 

with some C > 0, do > and a function g(p) > with (at least) g{p) = 0{p) for p — > 0, 
is well established (see Ivers &; James 1984) and is the prerequisite of the present section. 
We come back to this subject when discussing decay rates in section 5. 

The central tool of this section is the following higher-order-decay theorem for smooth 
solutions u of the system 

dtu — aAu + b • Vu + cu = f 
Au = 



u and Vu continuous 
u(x, •) — > 
u{- ,0) = u 



in G x ]R + . 

in G x R + , 
in IT x H + , 
for |x| — > oo, 
on G x {t 



0}. 



(3.2) 



Here, G is a bounded domain in H n (n > 3) with complement G = M n \G and (sufficiently 
smooth) boundary dG. Associated to (|3.2p is the eigenvalue problem 



—Au = Xu in G, 

Au = in G, 

u and Vu continuous in H", 

n(x) —7- for Ixl 



oo, 



(3.3) 



which is well-defined and has a lowest eigenvalue Ai that is positive (see appendix C). 

To measure the spatial smoothness of the coefficients a, b and c we use the notation 
a G C k (G x R + ) (k > 0), which means that all spatial derivatives up to order k are 
continuous and satisfy bounds of the forrr|§ \D k a\ := max^^ sup^ xR+ |Z) a a| < K for 

some K > 0. a G Cf(Gx]R + ) means that, additionally, dtu is a continuous function with 
corresponding bound. 

Theorem 1 (Higher-order-decay) Let a, b, c £ C k (G x H + ) (k > 0) with common 
bound \D k a\, \D k h\, \D k c\ < K and a > ao > 0. Let, furthermore, f satisfy the higher- 
order energy decay condition 



V / \D a f(^t)\ 2 dx<C r 



-2d,t 



t > 



(3.4) 



|a|<fc 



Here, a enumerates all derivatives D a of order < k. 
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with some constants Cf > and df > 0, and let u be a classical (i.e. u £ Cf(G x M+)) 
solution of \3.S\) . Then, the energy decay condition 

I \u{x,t)\ 2 dx <C e- 2dot , t>0 (3.5) 
Jg 

implies the higher- order energy decay 

^2 [ \D a u(x,t)\ 2 dx < Ce~ 2dt , t>0 (3.6) 

as well as decay in the maximum norm: 

^ sup|L> Q u(-,t)| <Ce~ dt , t>0, (3.7) 
\a\=l G 

where I < k + 1 — n/2. Here, d := minjdj — e,do,df}, where di := aoAi is the 'intrinsic' 
decay rate of \3. 2\) . and C and C are constants depending on G, the initial value u$, K , 
Co, Cf,n, k and e > 0. 

The proof of this rather technical result can be found in (Kaiser 2012b). Yet some remarks 
are in order: the existence of classical solutions of the required type depends on the regu- 
larity of the coefficients (as specified in the theorem), on the regularity of the initial value 
uq and on suitable 'compatibility conditions' (for details see (Kaiser & Uecker 2009)). So, 
as a rule, for sufficiently regular data, conditions (|3.4p and (|3.5p imply decay to arbitrarily 
large orders. The decay rate of u is bounded by the lowest of the three decay rates di, do 
and dj; the intrinsic one di, however, can only 'almost' be attained since typically C and 
C diverge (at most algebraically) in the limit e — > 0. 

The 5-dimensional formulation (|2,16p of the meridional problem is obviously of the 
type to which Theorem 1 can be applied. So, setting / = 0, n = 5 and k = 4 the 
remaining condition (|3.5p takes the form 

/ M(x,t) 2 dx = \S 3 \ [ (-^M(p,z,t)Yp 3 dpdz<Coe- 2dot , (3.8) 

Jg 5 Jg 2 V J 

where | S3 \ is the volume of the 3-dimensional sphere S3 . Condition (|3.8p is clearly implied 
by (|3.ip when observing the axis-condition of the radial function g(p). 
As consequences we have then (13 . 6|) : 

[ \D a M(x,t)\ 2 dx<Ce- 2dt , t>0 (3.9) 

|a|<5 Gs 

and (ET71) with / = 0, 1 and 2: 

5 

sup\M(-,t)\ +^sup|d :Ei .M(-,i)| < Ce~ dt , t>0, (3.10) 
G5 i = i G5 
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sup\d Xi d Xj M(;t)\ <Ce~ dt , t>0, (3.11) 

i,j=l G 5 

which yield uniform pointwise bounds on 



I B m | 



-d z Me p + -d p Me 2 
P P 



- \V(p 2 M)\ = \2MVp + P VM\ 
P 



and 

|J«| = |V x B m | = - |A«M| = p\AM\, 
P 

respectively. As M is harmonic in G5, continuous in H 5 and vanishes at infinity, Ai takes 
its maximum and minimum in G5. This conclusion holds likewise for each component of 
|x|V.A/f. In fact, |x|9 a:i A^ solves the elliptic equation with negative zeroth order term, 

A ~ 2 R2 • v - up)" = in ^5, 

and, thus, obeys a maximum principle, too. With (|3. 10|) we can, therefore, estimate for 
t > 0: 

5 

sup|B m ( Vj t)| < 2sup|M(-,t)| + J2sup\x\\d Xi M(-,t)\ 

H K 5 ~! R 5 

5 

< 2sup|A^(-,t)| + 5^Bup|x||9 a . i >t(-,*)| 

< max{2, R} C e~ dt =: C Bm eT dm \ 

where R is the radius of a ball enclosing the conductor. As J a is restricted to the conductor 
one obtains directly by (|3. 1 1 [) : 



(3.12) 



sup 

G 2 



-Ja(-,-,*) <sup|AM(-,t)| <C Ja e~ dmt , t>0. (3.13) 



According to Theorem 1 the decay rate bound d m is given by d m = min{dj — e, do}- in 
section 6 we identify the intrinsic decay rate di with the rate of free decay of the meridional 
scalar; thus, do < di. Comparing the decay of meridional field and meridional scalar we 
find, therefore, in the case do < di the same bound for both quantities and in the case 
do = di almost the same bound. Concerning the constants C*B m and Cj a we made no 
effort to determine the exact dependence on the various parameters listed in Theorem 1. 
So, numerical values could be very large. From an observational point of view, however, 
these can always be compensated by sufficiently large (compared to the decay time l/<i m ) 
periods of time. Finally, the high degree of regularity of the data (boundedness of fourth- 
order derivatives) required by Theorem 1 is due to the method of proof, in particular, the 
use of embedding results in H 5 neglecting the axisymmetry of the problem. Using more 
refined proof techniques the regularity requirements are expected to decrease. 
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4 Decay of azimuthal field and meridional current 



In this section the unconditional decay of the azimuthal field B a and of the meridional 
current J m is proved exploiting the divergence-character of the governing equation in 3 
dimensions as well as the nonsingular formulation in 5 dimensions. A central tool is a 
positive axisymmetric solution P of the following auxiliary problem 

d t P-V- (4 V (p 2p )) +V-(vP) = inGxR + , j 

n-VP = on9GxR + , > (4.1) 

P(. , 0) = P > on G x {t = 0}, J 

that differs from (|2.15p in a vanishing right-hand side in (|4.ip i and a Neumann boundary 
condition instead of a Dirichlet condition. All data G, rj, v and Pq are assumed to 
be axisymmetric, which implies that the solution is axisymmetric, too. The regularity 
requirements that ensure the existence of solutions are moderate; they are more detailed 
for the non-singular 5-dimensional formulation in appendix A. We mention here only 
that G G H 3 is a bounded domain with smooth boundary dG and exterior normal n, 
rj G C 1 (G x H_|_) and v € C(G x IR + ) with n • v\qq = 0. Moreover, we need the bounds 
|v| < K, |Vp • v/p| < K, \Vp ■ Vt]/p\ < K for some K > and r? > r/ > on G x M + . 

Theorem 2 (Positive auxiliary solution) Let P be a solution of the axisymmetric 
problem ft4-l\ ) with initial value Pq and let P_ , Pq be some positive numbers that bound 

Eq<Po< Po in G. (4.2) 
Then, positive bounds P and P exist such that 

P<P(,t)<P fort>0, (4.3) 
where P and P depend only on G, rjo, K and P and Pq. 

The proof of Theorem 2 consists in an adaptation to our situation of a powerful result by 
Lortz et al. (1984). The necessary changes are rather technical and are, therefore, deferred 
to appendix A. 

We compute next the time-derivative of the quantity f G A 2 /Pdx. To avoid singular 
coefficients we consider first a 'regularized' domain G e := {x G G : p(x) > e} (e > 0) 
bounded by the surface S e := {x G dG : p(x) > e} and the cylinder C t := {x G G : 
p(x) = e} such that dG e = S e U C e . In the limit e — > 0, C e shrinks to the line segment 
L := {x G G : p(x) = 0}. By (f2TT5]) i. (gZEjb and with the abbreviation / := B m • V(^/p) 
one obtains^] 

dx= {2 — dtA-— T d t P)dx 



dt.l Gc P J Ge \ P 1 P 2 



2 L * f dx + L { 2 ^ v ■ (? v(pl4) - - S v • (^ v ^ p ) - vp ) } d " 



3 For simplicity, dependence on t is suppressed in the subsequent calculations. 
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Integrating by parts in the last integral and using the boundary condition ([2.150 9 on S £ 
the velocity field v drops from the integral and we are left with 



p- 



T}P 



dx 



A 



A 



V d p A - —r d p P )ds - 2 / - — ds 



P 



p2 



r\ A 2 



c, 



< C P P 



In the limit e — > only the second surface integral survives with the result —An J L r\ A 2 IP d 
In summary we have in the limit e — > 0: 



2. 



f ^- dx = -2 f rjP 

Jg p Jg 



A 2 

dx — Air I r] — dz + 2 



dtJ G P 

By means of the bounds ([4.30 on P, the variational inequality 



G 



P 



fdx. 



(4.4) 



/ |V#| 2 dx + 27r / \g\ 2 dz > m / |g| 2 dx 
Jg Jl Jg 



for axisymmetric differentiable functions g vanishing on dG (see section 6), the bound 
([3.120 on B m and the bound |V(u^/p)| < if, ([4.40 can be estimated as follows: 



d_ f A 2 ^ 
dtJ G P 



dx < -2rj P 



< -ZvoPiP 



G 



/A\ 2 

A 



dx + 2irJ L ^ 2 dz}+2j G ^B m .v(J)dx 



1 



< -2r? /ii= 



G 

P /■ -4 2 



p ,dx + 2^C Bm e-^ pl/2 „ 



.4 



dx 



G 



P 



dx + 2KC Bm e" 



IGI 1 / 2 

pl/2 



P 



dx 



1/2 



/• ^ 2 K 2 Cl 
< -2(4 - e) / — dx + 5^ 



K ^bJ^I e -2d m t 



(4-5) 

In the third line we used Youngs's inequality to split the second term and introduced the 
azimuthal decay rate d a := rjfj,\P/P. Applying Gronwall's inequality on ([4.50 yields then 
the exponential decay of the quantity f G A 2 / P dx: 



G 



A(x,t) 
P(x,t) 



dx < 



G 



Po(x) dX + 



K 2 Cl \G\ e -2(d a -e)t _ e -2d m t 



2eP 



2(d m - d a + e) 



t > 0. 



Eliminating, finally, the auxiliary function P, we end up with the energy estimate 
„4 2 (x,t)dx < P [ dx < — / ^(x)dx e -2(*.-«)* 

^Bml^l — 2min{d a — e,d m }t (~t — 2 min{d a -e, d m }t 

+ P4e\d m -d a + e\ e - UA& 



G 



(4.6) 
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where Ca depends on the initial value Aq, the fraction P/P, K, C& m , \G\, d m — d a and 
e > (chosen such that d m — d a + e 7^ 0) . 

To obtain bounds on B a and J m = V x B a in the maximum norm we use once more 
Theorem 1, which holds literally also in the case that the 'boundary condition' ()3.2j) -? 3 4 
is replaced by a Dirichlet condition (Kaiser 2012b). Resolving the divergence form, the 
governing equation (|2. 17[) i of the azimuthal problem in H 5 takes the form of (|3.2[) i . So, 
Theorem 1 may be applied to solutions of (|2.17p . Choosing n = 5 and k = 3, the 
main conditions fj3.4|) and (|3.5p are satisfied by (j3.9j) with dj = d = d m and (|4.6p with 
do = mm{d a — e, d m }, respectively. Note that the intrinsic decay rate d{ of (|2.17p coincides 
with the azimuthal free decay rate (see section 6), thus d a < d%. By (|3.10p we obtain then 
the bound 

sup|.A(-,t)| +sup|VA(-,t)| < Ce- min{d ^ e ' dm}t , t > 
G s G 5 

and hence for all time 



sup 

G 



-B a (-,t) = sup\A(;t)\<C Ba e~ min ^~ e ' d ^ t (4.7) 

P Gs 



and 



sup|J m (-,t)| < 2 sup \A(-,t)\+R sup \VA(-,t)\ < Cjm e- min ^- € ' d ^ t . (A, 

G G5 G5 



The constants C^ a and Cj m depend, as in the meridional case, on bounds on the deriva- 
tives of the data rj and v up to fourth order; moreover on r]o, the conducting region G, 
the initial values Mq and Ao, and the 'balance parameter' e. All other dependencies can 
be eliminated by substitution; for instance the 'oscillation' P/P of the auxiliary func- 
tion P depends, according to Theorem 1, on no other parameters when the initial value 
Pq = const = 1 is chosen. 



5 Meridional decay rates 

In this section we present the example of a flow field, for which the meridional decay rate 
shrinks exponentially fast to zero with respect to the amplitude of the flow. For simpler 
flows, viz. purely radial or purely non-radial flows we prove the opposite behaviour, i.e. 
even for large amplitudes the meridional decay rates do not drop significantly below the 
corresponding free decay rates. These results are obtained by constructing suitable sub- 
and supersolutions for the dynamo problem. Their significance depends on the following 
maximum principle that is formulated according to our needs. 

Let G C H n be a bounded domain with smooth boundary dG and complement G. 
Let, furthermore, G be decomposed into a finite number of subdomains G\, . . . ,Gk by 
smooth hypersurfaces Ti, . . . , Ti such that G = (Ji=i Gi an d G% H Gj = for i ^ j. So, 
dGi may have parts in T := U?=i Tj and dG; it is smooth up to possible intersections of 
the bounding components Tj and dG. To take care of a possible symmetry axis S in the 
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(5.1) 



x n -direction we use the notation x = (p,x n ) with p = \p\ = ( Y^j=i x 'j) 1 ^ 2 - ^ denotes 
the elliptic operator L := aA + b • V + c with coefficients a, b, c satisfying the regularity 
requirements a, Va, ph, pc G f|Li <?°(Gi x [0, T]) n C°{G x [0, T]) and a G C°(G x [0, T]) 
with a > ao > 0. So, b and c are possibly unbounded at S and need not be continuous 
over T, whereas a is supposed to be continuous in all G. Functions u = u(p, z, t) then 
satisfy the system of inequalities 

dtu < Lu in Gi x (0, T) for % = 1, . . . , k, 

0<Lu in G x (0,T), 

u and V« continuous in H n x (0,T), 

w(p, -, •) < for p — > 0, 

w(p, ^, •) < for |x| — )• oo, 

n(-, -, 0) = no on G x {i = 0}. 

in a classical sense if all derivatives exist as continuous functions, i.e. u G n«=i x 
(0,T)) flC 2 (G x (0,T)). For such functions we have the following result 

Theorem 3 (Maximum principle) For functions u G (G x [0,T]) with T > 0, which 
satisfy 15. 1\) with c < in a classical sense, holds 

sup u < max u^. 

R"x(0,T] G 

Here, means the positive part of uq, i.e. (x) := max{uo(x), 0}, and u G C(G x [0, T]) 
implies that n takes continuously its initial value uq on G. A proof of Theorem 3 is given 
in Appendix B; we make here only a few comments: our result differs from similar ones 
in (Backus 1957) or (Ivers &; James 1984) in that several separate conducting regions are 
admitted. This allows us to consider discontinuous flows, which typically have no classical 
solutions in all G. Theorem 3 states, in case of a positive initial maximum, that u can 
never exceed this value and, in case of a negative initial maximum, that u can never exceed 
zero0 Theorem 3 holds likewise in the bounded case, where (|5j!])2,5 are replaced by the 
boundary condition n|aGx(o,T) ^ 0- 

In this section it is convenient to use polar coordinates (r, 9) in the meridional plane 
in the form p = rsin#, z = rcosd with r > and < 9 < ir. The symmetry axis is then 
represented by 9 = and ir, r > 0. Equation (12. 5p for the meridional scalar M = M(r, 9, t) 
takes then the form 

d t M - j] (df.M + ^ dg ( dgU) ) + v r d r M + -v e d e M = (5.2) 



sin V / / r 



4 We take the opportunity to correct the maximum principle in (Kaiser 2007), where we missed the case 
of a negative initial maximum: in Lemma 1 qo should be replaced by . All consequences we have drawn 
there from Lemma 1 remain untouched by this correction. 
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and the axis condition (|2.19p reads 



M(-,0,-) = O(sm 2 0) for 9^0 or vr. (5.3) 
Obviously, (r,9,<fi) constitute spherical coordinates in H 3 . 

5.1 Exemplary slow decay 

Let us consider the piecewise constant (discontinuous) velocity field v (in spherical coor- 
dinates) in Bi Cl 3 : 

f -3c 

Vr = \ 3c 

ve J — 9 c 
r sin \ 9 c 

with some constant c > 0; otherwise v is zero, in particular, v = in B i / S and = 0. 
Let M be a solution of the meridional problem (|2.14p with this flow field, with rj = 1 and 
with initial value Mq > 0. In general, M will not be a classical solution in G = B%, but M 
as well as — M will satisfy the premises of Theorem 3, especially (|5.ip with Si subdivided 
into -B1/3 and four spherical half-shells. Applying Theorem 3 on — M then yields 

inf M> mm{min{M (x),0}} = (5.5) 
R 3 x(0,T] xefii 

for any T > 0. 

In the following we construct a subsolution M_ of (15. 2p (i.e. the left-hand side in (I5.2p 
is non-positive) with decay rate d m on the half-annulus 1/3 < r < 1, < 6 < ir satisfying 
(15. 3p and the boundary condition M_\r=i/3 = M]r=i = 0. In IR 3 this region corresponds to 
the spherical shell SS := B\\Bi/^. The function M_— M then satisfies on SS the bounded 
version of (|5.ip . in particular, we have by (|5.5p . (M — M)|a<y? < 0. So, choosing M < Mo, 
Theorem 3 yields supgg y (p yi (M — M) < for any T > or, equivalently, M(-,t) > M_{-,t) 
for any f > 0, i.e. d m is an upper bound on the decay rate of M. 

To obtain a suitable subsolution M we make the ansatz 

M(r, 9, t) := f(r) g{9) e^™* (5.6) 

with non-negative functions / and g, and a constant d m > 0. The construction of / and 
g is based on the auxiliary function 

h(s) := e~ cs l 2 sinh^(c/2) 2 -ds, < s < 1, (5.7) 

where the parameters c and <i are related by 

V(c/2) 2 -d coth V(c/2) 2 -d = c/2, 



for 1/3 < r < 2/3 
for 2/3 < r < 1 

for < < vr/2 
for 7r/2 < < 7r 



o < e < vr, 



1/3 < r < 1, 



f > 



(5.4) 
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which defines a function d = D{c) that behaves asymptotically like 

d = D{c) ~ c 2 e~ c for c -> oo. (5.8) 
/i is the unique solution of the boundary value problem 

d 2 d d 

— h + c — h + dh = 0, h(0) = 0, —h(l) = 0. (5.9) 
ds^ as as 

Thus, with the substitution t := 1 — s, h satisfies the inequality 

d 2 d 
(1-t 2 ) — ^h-dl-t 2 )— h + dh>0, 0<t<l, 
at z at 

or, after the further substitution t := cos 9, 

s ' me ^(^-E^E9i) + csm9-^-g 1 +dg 1 > 0, < 9 < J, 
a9 \ sin 9 a9 J a9 2 

where g\{9) := h(t). At the boundaries we have <?i(0) = d/d9 gi(n/2) = 0. Introducing vg 
from (|5.4p and observing that d/d0 51 > we have, finally 

d/ld\ d , . „ „ „ 7T 1 



Sm ^l^d^ 5l J-^d^ 5l+ ^ 1 ^°' ° <0< 2' 3 <r<L (5 - 10) 

By reflection at = 7r/2 we find that g%{9) '■= 51 (tt — 9) satisfies inequality ()5.10|) in the 
region ir/2 < 9 < ir, 1/3 < r < 1. 

Concerning the radial part we start from (|5.9[) with d replaced by d := d m /9 — d and 
find by the substitution r := (1 + s)/3 for /i(r) := /i(s) the inequality 

r 2 ^h-r 2 v r ^h+r 2 d m h-dh > 0, 1/3 < r < 2/3 (5.11) 

and the boundary conditions /i(l/3) = d/dr/i(2/3) = 0. An analogous calculation yields 
/2(r) satisfying (|5.1ip on 2/3 < r < 1 with boundary conditions d/dr/ 2 (2/3) = / 2 (1) = 0. 

The functions fi^ and g\ : 2 constitute now the subsolution (|5.6p on their respective 
domains: (|5.10p and (|5.1ip imply (|5.2|) (as inequality) with 77 = 1 and one easily checks 
C 1 -smoothness over the subdomains, the asymptotic condition (|5.3p and the boundary 
conditions M(l/3, ■) = M(2/3, •>•) = 0. For large flow amplitudes c, the quantities cf 
and (i and hence d m obey (|5.8p . which manifests the (exponentially) slow decay of M. 



5.2 Purely radial and non-radial flows 

A super solution M (i.e. — M satisfies (]5.ip ) with decay rate d m that bounds a solution M 
at t = does so for all time, hence d m provides a lower bound on the decay rate of M . 
This follows by Theorem 3 applied on M — M and on —M — M. 
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We consider first purely radial flows, i.e. vq = 0, in an arbitrary conducting region G. 
Concerning v r the following bounds will play a role: 

-v r < 2 ci m/R for < r < R/2, ) 

0<e<Tr, t>0, } (5.12) 
v r < 2 c 2 7] /R for R/2<r <R, J 

where % is a lower bound on r\ and i? the radius of a ball enclosing G. For large flow 
amplitudes c\ , C2 a supersolution is then given by 

M r (r, 9, t) := f(r) sm 2 (6) e" d -*, (5.13) 

where d r m := 8r)o/R 2 and 

ai^e" ClX l3 /2 (6i2;) < x < 1/2, 

a 2 V^e C2:c (^3/2(622:) + 03/3/2(622;)) 1/2 < x < 1, 
x -1 x > 1 



with x := r/i? and bi := (c? — di) 1 / 2 (z = 1,2). and K v denote modified Bessel functions 
of order v and I' v and K' v their derivatives. Thus /(x) is a solution of 

x 2 f" -2{-l) i c i x 2 f + d lX 2 f -2/ = (* = 1,2) (5.14) 

on < x < 1/2 (i = 1) and 1/2 < x < 1 (i = 2), respectively. To ensure C^-smoothness 
we require /'(1/2-) = and /'(1/2+) = 0, which yield 

(1 - d)/ 3/2 (6 1 /2) + 61/3/2(61/2) = (5.15) 

and 

(1 + c 2 )(K 3/2 (b 2 /2) + 03/3/2(62/2)) + b 2 (K' 3/2 (b 2 /2) + 03/3/2(62/2)) = 0, (5.16) 

respectively. 03 is determined by the condition f'(l— ) = — /(l— ) and 01 and 02 are 
chosen such that / is continous. Relation (|5,15p determines a function d\ = D\{c\) that 
for large arguments decreases monotonically to the asymptotic value linic^oo D\{c\) = 80 
Similarly, ()5. 16j) determines a function D 2 with the same asymptotic property. Note that 
/ is monotonically increasing for < x < 1/2 and decreasing for 1/2 < x < 1, which 
allows us to use the bounds (|5.12p . Thus, M r satisfies (|5.2p (as inequality) in Br and 
hence in G. The remaining conditions in (|5.ip are easily checked. This verifies M r to be a 
supersolution with decay rate 8t]q/R 2 , which differs from the free decay rate ir 2 r]o/R 2 by 
less than 20%. This bound cannot be improved since for a piecewise constant flow field in 
a ball according to (|5.12p with c\/c 2 chosen such that d\ = d 2 , (|5. 13j) becomes an exact 
solution of the dynamo problem. 



5 This analysis has to be done with some care but it is elementary since I3/2 and -K3/2 can be expressed 
by elementary functions (cf. Abramowitz & Stegun 1972, p. 443). 
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For purely non-radial flows holds v r = 0. With the bounds on vq, 
—tvq < c smOrjo for < 9 < n/2, 

r < R, t>0, 

rvg < c s'mOrjo for ir/2 < 9 < 7T, 
a supersolution reads now: 

M nr (r,9,t) : =f(r)g(9)e- d - t , 

where := 7r 2 r/o/(4i? 2 ), f(r) := sin7rr/(2i?) for r < R and f(r) := 1 for r > R, and 
g{9) := h{\ - cos 9) for < 9 < vr/2 and g{9) := h{\ + cos 9) for vr/2 < 9 < vr with ft being 
the function (15. 7p . The verification of the conditions (15. ip is left to the reader. 



6 Incompressible fluids of constant conductivity 

When the flow field v is such that V • v = and the diffusivity r\ is constant, much 
larger decay rates can be obtained than without these restrictions. In fact, the decay 
rates of the meridional and the azimuthal field are those of axisymmetric free decay. If 
the conducting region is a ball, explicit calculations show that these decay rates coincide 
with the well-known poloidal and toroidal free decay rates, respectively. These numbers 
provide, moreover, lower decay rate bounds for any region enclosed by these balls. The 
basis of these results is a careful evaluation of the corresponding energy balances and of 
the associated variational problems. We start with the meridional case. 

Multiplying (|2.14p i by M, integrating over the regularized domain G t introduced in 
section 4, inegrating by parts and using the boundary condition n • v = on S e yields 



Id 

2 di 



/ M 2 dx = ( (r]AM - 2? ? — • VM - v • Vm] M dx 
Jc e Jg c v p j 

= -r] I \VM\ 2 dx + r] j M n • VM ds-rj j - ' Vp M 2 ds (6.1) 
Jg c Js e J s t P 

f f M 2 1 f 
-rj MVp-VMds + rj I ds + - Vp-vM 2 ds. 

An analogous calculation yields in G 6 := {x £ G : p(x) > e}: 

= 7] [ (AM — 2 — • VM] Mdx = —rj [ \VM\ 2 dx-r)[ Mn-VMds 
Jdc v P 1 Jg. Js. 

+r] [ n ' Vp M 2 ds - i] [ MVp-VMds + r] [ —ds, 
Js e P Jc e Jc c P 



(6.2) 



where we used the asymptotic conditions (|2.14p ^ and (|2.18p and the notation C e := {x G 
G : p(x) = e}. Summing up f)6. 1 1) and (|6.2p and letting e — > while observing (|2. 14|) -^ a 
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yields, finally, the energy balance 

— / M 2 dx = -2ri [ |VM| 2 drE. 
d* Jg Jr 3 

So, solving the variational problem 

/ R3 |VM| 2 dx 



inf 

0^MeS„ 



Ai 



f G M*dx (6 ' 3) 

where S m means the set of axisymmetric functions that satisfy ([2.140 3/ k. yields exponen- 
tial decay of the energy with rate 2dS^ n := 277A1 . The Euler-Lagrange equations associated 
to (j6.3p are obtained by a calculation analogous to that of (|6. 1 j) . (|6.2p . They amount to 
the eigenvalue problem 

-AM + 2 — • VM = AM in G, 

AM - 2 — ■ VM = in G, 

M and VM continuous in H 3 , 
M(x) -> , VM(x) = O(p) for p 0, 



P 



VM(x) -> 



for 



x 



00, 



(6.4) 



whose lowest eigenvalue Ai solves (|6.3[) (see Appendix C). Obviously, (|6.4p and hence Ai 
still depend on G but not on v; thus, Ai defines the free decay rate in G. 

When reformulated in H 5 , (16. 4p represents an eigenvalue problem for the Laplacian. 
This eigenvalue problem has explicitly been solved for the case that G C H n is a ball Br 
of radius R (see Appendix C). For n = 5 we find Ai = (tt/R) 2 , hence = ??vr 2 /i? 2 , which 
coincides with the poloidal free decay rate. 

Note that Backus considers in the ref. (Backus 1957) the more general variational prob- 
lem (28), which for p = k = 1 coincides with our problem (|6.3p . Neglecting axisymmetry, 
in particular, the axis-condition (]6.4[) ^i . Backus derives the Euler-Lagrange equations (30) 
and (31), which formally coincide with those describing poloidal free-decay modes in H 3 . 
Since in the axisymmetric situation the asymptotic condition at infinity is weaker and no 
zero- mean condition applies, he finds in balls Br the lowest eigenvalue (ir/2R) 2 , which sug- 
gests a decay time that is 4 times larger. This result, however, does not seem to be correct. 
In fact, the corresponding eigenfunction is spherically symmetric; so, its gradient does not 
vanish on the symmetry axis and, therefore, is not admissible in the representation (|2.2p 
of axisymmetric fields. 

Concerning the azimuthal case we can make use of the calculations already performed 
in section 4. With the notation introduced there the energy balance can be read off (I4.4p 
by setting P = 1: 

d 



dt 



G 



A z dx = -2r] / \VA\ z dx- Attt] / A z dz + 2 / Afdx. 



G 



(6.5) 



G 
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Neglecting the last term, the azimuthal variational problem is then 

L\VA\ 2 dx + 2ir [ T A 2 dz , x 

inf ^ ' r A9 , JL =:/ii, 6.6 

where S a means the set of axisymmetric differentiable functions vanishing on dG. The 
associated Euler-Lagrange equations read 

-V.(1V(^))=M inG, (6?) 
A = on dG, 

where the lowest eigenvalue hi solves the variational problem (|6.6p . This is again more 
transparent when reformulated in IR 5 . Problem (j6.7H becomes then the well-known eigen- 
value problem for the Laplacian with Dirichlet boundary condition: 



-AA = nA in G 5 , 
A = on dG 5 . 



(6.8) 



In balls Br the eigenvalues are explicitly known; in particular, there is [i\ = (i\/R) 2 , 
where i\ is the lowest positive zero of the first spherical Bessel function j\ (cf. Folland 
1995, p. 106). 

For the azimuthal field energy, the decay estimate (|4.6p holds then with P = Pq = 
P = P = 1 and azimuthal decay rate d® := rjfMi. In balls we have d® = r](i\/R) 2 , which 
coincides with the toroidal free decay rateH Note that (|4.6p contains the possibility that 
in the presence of a meridional field, the azimuthal decay is also dominated by the (lower) 
meridional decay rate. 

Pointwise decay can again be obtained by means of Theorem 1. In the azimuthal case 
without meridional field pointwise decay holds with almost the free decay rate since 
the five-dimensional energy (which is required by Theorem 1) is dominated by the three- 
dimensional energy. This is not so in the meridional case. Only when using additional 
information, five-dimensional energy decay (with a lowered decay rate) can be established: 
splitting the conducting region G into G e and its complement and using (13. ID we can 
estimate 

- \M f (±*,f- ^ 2 ~2, 

27T 7 G Vp2- 



M(x,tydx =h^r I (— Af(x,t)) p z dx 

G 5 

S 3 \ 1 



/ M(x,t) 2 dx + C 2 [ g( P ) 2 \dx 

JG e JG\G e P 



~ 2tt e 2 
" 2e 2 + 2 



6 This confirms again the assertion of Braginskii (1965), although he seems not to have considered the 
correct variational problem. Backus (1957) finds with ingenious but not completely rigorous arguments 
again a lower decay rate, viz. r](n/R) 2 . 
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Setting e := exp(— we thus have 

f M(x,t) 2 dx < Ce" 4 ' 

and hence by (|3.12|) pointwise decay of the meridional field with (at least) half the free 
decay rate (which is supposedly not optimal). 

Finally, one should note that eigenvalues of (|6.8p are monotonous with respect to the 
domain, i.e. G C G implies fj, > Jl, where \i and Jl belong to G and G, respectively (cf. 
Courant & Hilbert 1961, p. 409). The eigenvalues A of (|6.4p share this property (see 
Appendix C). So, the decay rates df^ and d® in arbitrary (axisymmetric) domains G can 
always be estimated by those for balls enclosing G. 



Appendix A 

This appendix expounds the necessary changes in the proof of (Lortz et al. 1984, Theorem 
2) to be applicable to Theorem 2. First we rewrite system (|4.ip in H 5 . In view of (|2.17p 
one obtains 

d t P -V ■ (rjVP) + V ■ (v$P) + cP = mG 5 x (0,T), 

n-VP = on 3G 5 x (0,T), } (A.l) 

P(. , 0) = P > on G 5 x {t = 0} 



with c := —2(d p 7]/p + v p /p) and arbitrary T > 0. System (jA.ip differs from system (|2.2p 
m (Lortz et al. 1984)0 by the zeroth-order term cP. Classical solutions are guaranteed 
for both systems by the sufficient conditions dG G C 3 , rj and G C^(Gs x [0, T}) with 

n • v^laoxflVTI = 0, c G Cj(G 5 x [0,T]) and P G C 3 (G 5 ) with n • VP \ dG = 0, which 
imply the precise Holder conditions formulated, e.g., in (Ladyzenskaja et al. 1968, chap. 
IV, Theorem 5.3). The proof of (Lortz et al. 1984, Theorem 2) refers to such classical 
solutions. 

The governing equation (LMS2.2a) enters the proof two times, viz. in the derivations of 
the inequalities (LMS3.1) and (LMS3.12). We show in the following that the zeroth-order 
term in (jA.ip i does not invalidate these inequalities, it merely modifies the bounds that 
depend then, additionally, on c. 

Modified proof of (LMS3.1): Multiplying (]A~l|) i by P^ 1 with 7 ^ 0,1, integrating 

(5) 

over G, integrating by parts and using the boundary conditions for P and Vm yields 

[ P^dx = [ d t PP 1 dx= ! {V- CnVP- v^ ) P)P 7 ~ 1 -cP 7 )dx 
7 at J G J g J g ^ a 2 ^ 

{(1 - 7 )(r ? P^ 2 |VP| 2 - P 7_1 v^ • VP) - cP^}dx. 



7 Equation numbers referring to (Lortz et al. 1984) are henceforth marked by the suffix 'LMS', thus, 
'(LMS2.2)'. 
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With the bounds r\ > 770, |v^| < K and |c| < K, and by Young's inequality, (|A.2[) can be 
estimated as follows: 



dx 



—r - — t-4- / ^ 7 da; < -r/o / P 7 ~ 2 |VP| 2 dx + if / P 7 ~ 1 |VP|dx + , - , / P" 
7(7 - 1) 7 G 7 G y G I7 - 1| Jg 

<_% f P ^\ VP fdx +(— + -?—) [ P 7 dx. 
4 Jg v % |7-1k7g 

Finally, multiplying by 4>(t), integrating over [t±, £2] and integrating by parts yields (LMS3.1) 
with the constant K 2 jr\Q replaced by K 2 /770 + K/Y) — 1|. In the course of proof the pa- 
rameter 7 takes infinitely many values; however, always holds I7 — < n + 1 if n > 3, 
in particular, (7 - l]" 1 < 6 in R 5 (Lortz et al. 1984, p. 688 and p. 690). Thus, (LMS3.1) 
holds with an enlarged constant. 

Modified proof of (LMS3.12): Multiplying (jXTTl i by P" 1 and integrating over G yields 
analogously: 

^ / logPdx = / {V-(??VP-v^ ) P)p- 1 -c}dx 
dt Jg Jg 

= / ?7|VlogP| 2 dx - / V(logP) • dx - / cdx 
Jg Jg Jg 

>lva [ |VlogP| 2 dx-i — \G\-K\G\, 
2 J G 2 Vo 

which is (LMS3.12) with the constant K 2 \G\/(2r]o) enlarged by K\G\. 



Appendix B 

Proof of Theorem 3: Let M denote sup R ™ x (o,T] u i which is finite because of (|5.1|k . The 
case M < is trivial: by continuity we have sup^uo — an d hence M < = maxgjij. 

To prove the case M > we make use of the elliptic maximum principle and the elliptic 
boundary derivative theorem for classical solutions as formulated in (Protter & Weinberger 
1984, Theorems 6-8 at p. 64ff) and the parabolic maximum principle for weak solutions 
as formulated in (Lieberman 1996, Theorem 6.25 and Corollary 6.26 at p. 128). We show 
first that supgu(-,t) < max^ «(•,£) for any t G (0, T] and, second, that max^ x j or jU is 
taken at t = 0, which proves our assertion. 

1) Let to ^ (0>^1 an d M to := supgu(-,io)- Because of f)5. s we find a ball B with 
complement B and a solid cylinder SC around the symmetry axis with complement SC 
such that u(-,to) < Mf in B U SC. Applying the elliptic maximum principle on the 
bounded and regularized domain GflBfl SC we find Mt be attained at some point 
(xo,to) £ 9G = dG. Furthermore, as u 7^ const in G and dG being smooth the boundary 
derivative theorem is applicable with the result 

n-Vu| (x0jto) <0, (B.l) 



24 



where n denotes the exterior normal with respect to G at (xo,io)- By (|5,ip ^ condition 
(jB.ip implies w(x, to) > M to for some x G G and hence M > M to for any to, i.e. 

M > u\ dGx ^ T] . (B.2) 

2) As u need not be a classical solution in all G x (0, T) a parabolic maximum principle 
for weak solutions is better suited to complete te proof. In order to check that u satisfies 
(|5.ip i in the weak sense of (Lieberman 1996, p. 100 above) we write L in the form 

L = V-(aV) + (b-Va)-V + c (B.3) 

and choose again a solid cylinder SC with boundary C such that u < M in SC n G, which 
implies, in particular, 

M > u|(cnG)x(o,T] • ( B - 4 ) 

Let v G C 1 ((5C nG)x [0, T]) be a non-negative testfunction that vanishes at d(SCf~) G) x 
[0, T]. Integrating vLu over (SC n Gj) x (0,T), integrating the first term by parts and 
summing over % = 1, . . . , k reveals that u satisfies indeed (I5.ip i in (SC flG) x (0, T) in the 
weak sense. Note that the regularity assumptions on u and the coefficients are such that 
the boundary terms on T cancel each other; otherwise, boundedness of the coefficients 
m (IB~3l) and the sign condition on c are the only further prerequisites of the parabolic 
maximum principle. Therefore, M is taken at the 'parabolic boundary', which means by 
(fET2|) and ([FTij) that M is taken at t = 0. 

The bounded version of Theorem 3 is obviously proved by 2) and (|B.2p . which holds 
now by virtue of the boundary condition u\qc X (o,t] ^ 0- 



Appendix C 



In this appendix we consider the (not necessarily axisymmetric) eigenvalue problem in IR" 

(n > 3) 

— Au = A u in G, 



Au = in G, 

u and Vu continuous in JR n , 
u(x) —7- for Ixl 



oo. 



(C.l) 



We give a variational characterization of the eigenvalues, demonstrate their monotonicity 
with respect to the bounded domain G C H n and compute them explicitly for balls Br. 
The axisymmetric case n = 5 is equivalent to the meridional eigenvalue problem (|6.4p . 

For sufficiently smooth boundary the system (IC.ip is known to have a countable set 
of (classical) eigensolutions {(iij,Aj) : i G IN} with positive real eigenvalues Aj, which 
are monotonically increasing with increasing i. The eigenfunctions constitute a complete 
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orthonormal set in L 2 (G), the set of square-integrable functions on G (cf. Kaiser & Uecker 
2009). 

Concerning the variational character of the eigenvalues observe that (jC.ip are just the 
Euler-Lagrange equations of the variational problem 



inf m ' ' =:Ai 



(C.2) 



with u varying in the set S" of functions satisfying (|C1|1 ^ a. So, Ai is the lowest of the 
eigenvalues of (jC.ljl and is realized by the associated eigenfunction u\. Neglecting ques- 
tions of degeneracy the next larger eigenvalue A2 then solves (|C.2|) if S is restricted to the 
orthogonal complement of u\. Repeating this procedure any eigenvalue has a variational 
characterization of type (|C.2p . 

The monotonicity of the eigenvalues with respect to G is now a direct consequence of 
(IC.2j) . since the denominator is always increasing with increasing G. 

If G = Br C H n , eigenfunctions and eigenvalues can explicitly be expressed in terms 
of spherical Bessel functions j u and n-dimensional spherical harmonics Y Vfl (cf. Folland 
1995, p. 98ff): 




/ -u+k— 1\ 



U/J, 



n) 

x / 



2—n—u 



Me 



in Br, 
in Br, 



} (C3) 



v+k-1 



/Rf 



where k := (n — 3)/2, I G W, v € INo, // G {1, . . . , -Dt/} and is the Zth positive zero of j u . 
()C.3P generalizes the case n = 3 (Kaiser & Uecker 2009, Appendix D) to arbitrary n > 3. 
That the set {u^} is an orthonormal set of eigenfunctions of (|C.lj) may be checked by 
explicit calculations. The completeness follows by analogous considerations as in (Kaiser 
& Uecker 2009, Appendix D). 

For the meridional problem (|6.4[) the lowest eigenvalue is by (jC.3p ? A^qq = (i®/R) 2 = 



(ir/R) 2 . The corresponding eigenfunction u^j is by (|C.3P i spherically symmetric , hence 



P 2u wo ^ s axisymmetric and satisfies the axis condition ([67 
in the meridional problem. 



1, 



and is, therefore, admissible 
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